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In this talk presented at the 5th Rencontres du Vietnam 2004, we discuss our preliminary 
investigations into voids of primordial origin. We show that if voids in the cold dark matter 
distribution existed at the epoch of decoupling, they could contribute significantly to the 
apparent rise in cosmic microwave background (CMB) power on small scales detected by the 
Cosmic Background Imager (CBI) Deep Field. Here we present the preliminary results of 
our improved method for predicting the effects of primordial voids on the CMB in which we 
treat the voids as an external source in the cold dark matter (CDM) distribution, employing a 
Boltzmann solver. Our improved predictions include the effects of a cosmological constant (A) 
and acoustic oscillations generated by voids at early times. We find that models with relatively 
large voids on the last scattering surface predict too much CMB power in an Einstein-de Sitter 
background cosmology but could be consistent with the current CMB observations in a ACDM 
universe. 



1 Introduction 

Analyses of surveys such as the 2-degree Field Galaxy Redshift Survey and the Sloan Digital 
Sky Survey indicate large volumes of relatively empty space, or voids, in the distribution of 
galaxies . In the standard hierarchical model of structure formation, gravitational clustering is 
responsible for emptying these voids of mass and galaxie^. However, standard cold dark matter 
(CDM) model simulations predict significant clumps of matter within voids that are capable of 
developing into observable bound objects and are not seen in the surveys^^^^. Peebles gives an 
in-depth discussion of the contradictions of this prediction with observational. He argues that 
the inability of the CDM models to produce the observed voids constitutes a true crisis for these 



models. Additionally, deep field observations from the Cosmic Background Imager^ (CBI) show 
excess power on small angular scales, I > 2000, in the cosmic microwave background (CMB). 

It may be possible to explain these 2 observations by postulating the presence of a void net- 
work originating from primordial bubbles of true vacuum that nucleated during inflation^ In 
this scenario, the first bubbles to nucleate are stretched by the remaining inflation to cosmologi- 
cal scales. The largest voids may have had insufficient time to thermalise before decoupling and 
may persist to the prese nt d av. Such primordial voids are predicted to produce a measurable 
contribution to the CMbHH. 

In this paper, we discuss the results of Griffiths et al. (2003)^^ in which we develop a general 
method to approximate the void contribution to the CMB that allows the creation of maps and 
enables us to consider an arbitrary distribution of void sizes. We show that if the voids that we 
see in galaxy surveys today existed at the epoch of decoupling, they could contribute significant 
additional power to the CMB angular power spectrum between 2000 < £ < 3000. Further to this 
work, we describe how we improve our predictions. We include the effects of the cosmological 
constant as well as oscillations in the matter-radiation fluid that may be generated by primordial 
voids on scales up to the sound horizon. For full det ails of our improved methodology refer to 
our recently submitted manuscript on the subject^l. 



2 First Approximation 

2. 1 Void parameters 

We model the voids seen today as spherical underdensities of 5p/ p = —1. Each void is bounded 
by a thin wall containing the matter that is swept up during the void expansion. This forms a 
compensated void. This means that the overall cosmology is that of the background universe, 
since a compensated void does not distort space-time outside of itself. This is an extremely 
important property, as it allows us to place many voids into a universe without the worry that 
they might influence each other. As a first approximation we take the bac kground univ erse to 
be an Einstein-de Sitter (EdS) cosmology. Conservation of momentum^^ or energy^l can be 
used to show that compensated voids in an EdS background cosmology will increase in radius 
r„ between the onset of the gravitational collapse of matter at equality and the present day such 
that 

rj,(r/) oc r/^ , (1) 

with (3 ~ 0.39 and where r] is conformal time. ,— ,— , 

Motivated by the extended inflationary scenario^^^, we assume a power-law distribution 
of bubble sizes greater than a given radius r of the form, 

iVB(> r) ocr"". (2) 

Tvpic ally, extended inflation is implemented within the framework of a Jordan-Brans-Dicke 
theory'-'^. In this case, the exponent a is directly related to the gravitational coupling uj of the 
scalar field that drives inflation, 

" = =' + ;7Tv5' 

Values of w > 3500 are required by solar system experiments'^. We take the limit of large w, 
leading to a spectrum of void sizes with a = 3. 

It has been shown that a distribution of void sizes that predicts the existence of arbitr arily 
large inflationary voids will cause significant effects on the CMB that contradict observation^^. 
We assume that the mechanism creating the voids imposes an upper cut-off on the size distribu- 
tion. A possible mechanism for this cut-off could be that the tunneling probability of inflationary 
bubbles is modulated through the coupling to another field. 




Figure 1: Top: The CMB anisotropies produced by the fiducial EdS void model (solid line) compared to the 
primary ACDM CMB anisotropies (dashed-triple-dotted). Also plotted are the sum of primary and void contri- 
butions (dotted) as well as the fluctuations induced purely by voids on the last scattering surface (dashed) and 
by those between last scattering and today (dashed-dotted). We show the "standard" cosmological concordance 
model: of course a combined analysis of primary and void-induced fluctuations would select a different cosmology 
for the primary contribution. Bottom: Example models depicting a range of void contributions to the CMB 
fluctuations. The models plotted are q = 3, rmax = 25 Mpc//i and Fv = 0.4 (solid line), a = 3, Tmax = 40 Mpc/h 
and Fv = 0.4 (dotted), a = 3, rmax = 40 Mpc/h and F^ = 0.2 (dashed) and a = 6, r^ax = 40 Mpc/h and 

K = 0.4 (dashed-dotted). 



As well as avoiding the well known problems associated with arbitrarily large voids existing 
at the epoch of decoupling, this assumption allows us to match the observed upper limit on void 
sizes f rom the galaxy redshift surveys. We choose the a,vGra,gG Vcl1ii6 that is found., ''^max — 

25 

Mpc/Zi'"''''^'. The minimal present void size is also chosen to agree with redshift surveys, rmin = 10 
Mpc/h. 

We normalise the distribution by requiring that the total number of voids satisfies the ob- 
served fraction of the universe filled with voids today, F^j. Redshift surveys indicate that ap- 
proximately 40% of the fractional volume of the universe is in the form of voids of underdensity 
5p/p < —0.9^^, ie. Fy ~ 0.4. The positions of the voids are then assigned randomly, mak- 
ing sure that they do not overlap. In order to speed up this process, we consider only a 10° 
cone. This limits our analysis to i > 100, which is satisfactory for our purpose since the main 
contribution from voids is on much smaller scales. 



2.2 Stepping through the void network 

Refer to Griffiths et al. (2003) 'H- for a more detailed description of our methodology. We ray 
trace photon paths from us to the last scattering surface (LSS) for the 10° cone in steps of 1'. 
Each void in the present day distribution that is intersected by the photon path is evolved back 
in time according to equation to determine whether the photon encounters the void. 

If a p hoton intersects a void between us and the LSS, we compute the Rees-Sciama (RS) 
effect due to the deviation in the redshift of the photon as it passes through the expanding 
void and the lensing effect due to the deviatio n in its path. If a photon intersects a void on the 
LSS, we calculate the Sachs-Wolfe (SW) effect^^ due to the photon originating from within the 
underdensity. We take into account the finite thickness of the LSS, which suppresses the SW 
effect for small voids, by averaging the contribution from a number of photons originating from 
a LSS of mean redshift 1100 and standard deviation in redshift 80. We also calculate the partial 
RS effect that arises due to the expansion of the void on the LSS as the photon leaves it. 

Once the photon has reached the last scattering surface, we know the variation of its tem- 
perature as well as its p osition on the LSS and can create a temperature map. We then use 
a flat sky approximationl^^^ to obtain the spectrum of the anisotropies (see Fig. P). We 
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Figure 2: The predicted CMB power from a uniform distribution of equally sized voids in EdS (dotted) and 
ACDM (solid) background universes. The voids have a size of 40 Mpc/h and a filling fraction of 40 %. The void 
contribution is computed by ray tracing a void network. The effect of A is to suppress the effect of the voids 
and move the peak in the power spectrum to smaller angular scales since they will appear smaller on the last 

scattering surface. 

point out that primordial void parameters are still poorly constrained by both observation and 
theory. The bottom panel of Fig. ^ shows a few further example models. 

For a power-law size distribution (as motivated by the inflationary scenario), large voids 
become rarer as a is increased. Therefore, since void analyses of redshift surveys only sample 
a fraction of the volume of the universe, there may exist voids of larger rmax than currently 
observed. Models with high rmax voids in an EdS background cosmology tend to predict too 
much power on scales i ~ 1000. However, if we take inflationary models with a > 6 then the 
peak moves to larger i and the total power drops. The filling fraction mainly adjusts the overall 
power. 

3 Improved Prediction 

3.1 Void Evolution in a ACDM background cosmology 

The addition of a cosmological constant (A) to the cosmology is expected to slow down the 
conformal evolution of the voids at late times. We test this hypothesis by simulating the comov- 
ing evolution of a single void of radius 25 Mpc//i starting at z = 1000 in a (100 Mpc/hf box 
containing 64^ particles. We find that the void size evolution in a ACDM background deviates 
from that of the EdS scaling solution at late times as expected. However, the final radius is 
underestimated by less than 2%. The EdS void scaling relationship given by equation with 
/? = 0.39 is therefore a good approximation for a void in a ACDM background. 

Since a ACDM universe evolves for substantially longer, we would expect the voids that we 
see today to appear smaller on the last scattering surface for this cosmology than in our first 
approximation EdS background. We would therefore predict that voids in a ACDM universe 
will have a more suppressed effect on the CMB than we have first approximated. This can 
be modelled by taking the horizon size of our cone of voids to be that of a ACDM cosmology. 
Fig. 121 compares the predicted power from a uniform distribution of equally sized voids in EdS 
and ACDM background cosmologies. The overall contribution to the CMB from the voids is 
suppressed as expected and the peak is moved to smaller angular scales since the voids now 
appear smaller on the last scattering surface. 

3.2 Modelling acoustic oscillations 

So far we have studied voids which do not interact with their surroundings. Specifically, they 
do not perturb the matter-radiation fluid themselves. While this is a reasonable approximation 
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Figure 3: A "zoom" onto the sound horizon showing the behaviour of the perturbations in real space at last 
scattering (LS) from a single void. The void, which has a radius of 40 Mpc today, has a size of 9 Mpc at 
LS. The black solid curve is the input perturbation "l>s(a;,7;Ls)/2, while the red dotted curve is the total metric 
perturbation $/2. The blue small dashed curve shows —2^ = 2$. The cyan large dashed curve is the temperature 
perturbation from the photons, Dg/A. The baryon velocity is shown as the magenta dot-dashed curve. The green 
small dashed-large dashed curve is Dg/A ~ 2"!?, which is very similar to $3/2. 

for the RS effect after recombination, it has been pointed oufl^that this is not true before last 
scattering. A void could therefore set up oscillations in the matter-radiation fluid on scales up 
to the sound horizon which might be clearly visible in the power spectrum. 

The standard way to include the full behaviour of matter and radiation perturbations is to 
solve the Boltzmann equation numerically. But it is not possible to model a void self-consistently 
inside the Boltzmann solver. A similar problem was encountered a few years ago while inves- 
tigating topological defects . As opposed to defects, which are entirely external sources, the 
voids represent a perturbation of the dark matter itself. We have chosen to model this by 
making the approximation that the dark matter is completely decoupled from the rest of the 
universe. Therefore, the cold dark matter acts only as an external seed where the perturbations 
are concerned (of course it is taken into account for the background quantities). We insert our 
source into a modified version of CMBEasyl^ and then write out a snapshot of the fluids at the 
time of last scattering and Fourier transform them back. 

Fig. 121 shows a close up of the perturbations for our model of void formation around the 
sound horizon. We see that the amplitude of the sound waves is only about 1% of the amplitude 
of the temperature perturbations inside the void. For smaller voids, the relative amplitude will 
be even less. This seems quite unimportant, however, the void itself covers a surface of only 
about 9^7r (Mpc//i)^ on the LSS, the sound horizon in our flat matter dominated universe of the 
example case is at lllMpc//i and the sound waves extend out to about 100Mpc//i, ten times 
further than the size of the void. Therefore, the total power in the fluctuations is approximately 
100 times larger than expected from their amplitude and is comparable to the the power from 
the void predicted by our ray tracing method. Indeed, we see oscillations appear on large angular 
scales in the angula r po wer spectrum (see Fig. EJ. For further details of our method, please refer 
to Ord et al. (2005)Il2l. 

4 Summary 

Fig. 0] compares the void Ci from ray tracing and the Boltzmann approach for a uniform dis- 
tribution of 40 Mpc//i sized voids. The overall effect from the first approximation ray tracing 
method appears to be suppressed using the improved Boltzmann solution. Evidence of acoustic 
oscillations generated in the photon-baryon fluid can also be seen on large angular scales. The 
RS effect due to alone, using the same approximations in both calculations, is in agreement 
for both methods as expected. 
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Figure 4: We compare the Ce from ray tracing and the Boltzmann approach. The dotted curve shows the total 
ST for a ACDM model with voids of a size of 40 Mpc/h and a filling fraction of 40 %, computed by ray tracing a 
void network. The solid curve shows the same but computed with a Boltzmann code. The small and large dashed 
lines show the RS effect due to alone for both methods and with the same approximations, these two curves 

should coincide. 
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Figure 5: Top: The CMB anisotropics produced by a ACDM void model (solid line: a = 3, rmax ~ 55 Mpc/h and 
Fu = 0.4) compared to the primary ACDM anisotropics (dotted) . Also plotted are the sum of primary and void 
contributions (dashed). We show the "standard" cosmological concordance model. Again a combined analysis of 
primary and void-induced fluctuations would select a different cosmology for the primary contribution. Bottom: 
The large dashed curve shows the fiducial EdS model of Fig. Q (a = 3, rmax = 25 Mpc//i and — 0.4) using 
our ray tracing method. The same model is shown for our Boltzmann solution with EdS (small dashed) and 
ACDM (dotted) background universes. The solid bold line is the same ACDM void model for rmax ~ 55 Mpc/h 

maximally sized voids as in the panel above. 



The suppression that is evident using the Boltzmann solution over the ray tracing approach 
is due to a number of factors. In the ray tracing method we assume that the density contrast 
5 = —1 and we take {a' /a)'^ to be Close to radiation domination (ie. at early times, just 

after last scattering, when we get the biggest effect) both these approximations act to increase 
the ray tracing result with respect to the more accurate Boltzmann prediction, together by 
more than a factor of 2. The ray tracing method also neglects the Doppler contribution from 
the baryon velocity, which suppresses the result further. Finally, the finite thickness of the 
LSS and Silk damping were not fully taken into account in our first approximation ray tracing 
method. 

As discussed in subsection 12.21 there may exist voids of larg er TjYiax than currently observed. 
Models with a = 3 and high rmax voids in an EdS background cosmology tend to predict too 
much power on scales i ~ 1000. Our results show that this is less of a problem for a ACDM 
universe. Furthermore, since the Boltzmann solution predicts the CMB power from voids to be 
lower than implied by the ray tracing estimation, relatively large voids on the last scattering 
surface may be consistent with the current CMB data, even for void distributions with low 



values of a (see Fig. [SJ . 

Experiments such as CBI are able to directly probe small angular scales and constrain void 
parameters. We will present a Markov Chain Monte Carlo constraints analysis of a wide range 
of void models in a future paper. We will further constrain models that are compatible with 
CMB observations using cluster evolutionl^^ and also investigate the non-Gaussian signal of void 
models that are compatible with the observations. 
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